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Abstract 

We prove the existence of scattering solutions for multidimensional 
magnetic Schrodinger equation which belong to the weighted Sobolev 
space H 1 _ s (R n )(n = 2, 3) with some 5 > \. As a consequence of this we 
formulate the direct Born approximation for the magnetic Schrodinger 
operator. Possible connections with inverse problems (inverse scatter- 
ing Born approximation) are discussed. 

1 Introduction 

The main goal of present article is to justify the application of the classical 
direct scattering Born approximation for the magnetic Schrodinger operator. 
The direct Born approximation is known as the most applicable approximate 
method in the numerous practical problems. It is also known that the inverse 
scattering Born approximation is well-defined and perfectly works (as the 
mathematical tool) in the case of linear and nonlinear Schrodinger operators 
and for all types of scattering data: full scattering, backscattering, fixed 
angle scattering and fixed energy scattering. For some scattering data it is 
possible to get the uniqueness and reconstruction procedure while for some 
data we are able to reconstruct singularities and jumps of unknowns even 
when there is no uniqueness. We mention here the results of Paivarinta and 
Somersalo [10] , Nachman (7] , [8] , Sun and Uhlmann (26] , Isakov and Sylvester 
[I], Paivarinta, Serov and Somersalo [To], Paivarinta and Serov [TT], [T3], [2], 
Ola, Paivarinta and Serov [9], Ruiz [TS], Ruiz and Vargas [19J, Paivarinta 
and Serov [12], Reyes [IB] . Serov [20], Serov and Harju [21], [22], Serov and 
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Sandhu [23J, Lechleiter [5], Reyes and Ruiz [T7] and some others. The main 
point of all these results is the precise calculation of the first (quadratic) 
nonlinear term in the Born series. For the magnetic Schrodinger operator 
the direct scattering problem (i.e., existence of the scattering solutions) as 
well as the inverse scattering Born approximation are not familiar at all. The 
big interest to this problem is connected to the fact that the knowledge of the 
scattering amplitude with backscattering data allows us to obtain essential 
information about the unknowns. 

We consider the magnetic Schrodinger operator 

H = -(V + iW(x)) 2 + V(x)-, xeR n , (1.1) 

in dimensions n = 2,3, where the coefficients W(x) and V(x) are assumed 
to be real- valued. We assume generally that W(x) G Lf(R n ) and 

VW(x) E L p s (R n ), n = 3, 3 < p < oo; n = 2, 2 < p < oo (1.2) 

and 

V(x) G L p s (R n ), n = 3, 3<p<oo; n = 2, 2 < p < oo, (1.3) 

where 5 > — -. Here L p a denotes usual weighted Lebesgue space and 
Sobolev space Wp a is understood so that / belongs to Wp^{R n ) if and only 
if / and V/ belong to L p a (R n ). For the case p = 2 instead of the symbol W\ a 
we use the symbol H\. 

It is well-known that under these conditions for the coefficients of the 
magnetic Schrodinger operator the following Garding's inequality holds: 

(Hu,u) L 2 {Rn) > v\\Vu\\ 2 L 2( Rn) - C\\u\\ 2 L2{Rn) , 

where < u < 1,C > 0. This inequality allows us to define symmetric 
operator H by the method of quadratic forms. H has a self-adjoint Friedrichs 
extension with the domain (in general) 

D(H) = {/(X) G Wl{R n ) : Hf(x) G L 2 (R n )}. 

In our particular case it is possible to prove that actually 

D(H) = W 2 (R n ). 

In the scattering theory the main role are played by the special solutions 
of the equation 

Hu(x) = k 2 u(x) 
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which are of the form 

u{x) = Uo(x) + u sc (x), 

where uo(x) = e lh ^ x ' 6 ' is incident wave with direction 9 E S™ 1 ' 1 and the 
scattered wave u sc (x) satisfies the Sommerfeld radiation condition at the 
infinity, i.e. 

n-i / du sc (x) \ 
lim r 2 I — _ iku sc (x) =0, r = \x\. (1.4) 



r^+oo \ dr 

In this case the total field u satisfies the so-called Lippmann-Schwinger equa- 
tion 



u = u + J G+{\x - y\) (iV{W(y)u) + iW(y)Vu - q(y)u") dy, 



;i.5) 



B" 



where q = \W\ 2 + V and G^ is the kernel of integral operator (— A — k 2 — z0) _1 . 
Using the representation u = u + u sc we rewrite this integral equation (1.5) 
only for scattered field 11 sc clS 

u sc = u + J G+{\x - y\) (iV{W{y)u sc ) + iW(y)Vu sc - q{y)u s ^ dy, (1.6) 
where u is equal to 



u {x) = / G+(\x - y\) (iV(W(y)u ) + iW(y)Vu - q(y)u ) dy. (1.7) 



We use the following results of Agmon [2] (see Remark 2, Appendix A): 
1 



^ „., HHlsiRn) + + \k\\\f\\Ll_ s (R«) < 4(A + k 2 )f\\ Lj{Rn) , \k\ > 1, 

where 5 > \ and H 2 s (R n ) denotes the weighted Sobolev space. As a conse- 
quence we have (for fixed k) that 

ii(-a - k 2 - xorviifl^) < c(k)\\f\\ Lj{Rn) , 

and uniformly in |fc| > 1 we have that 

\\(-A-k 2 - i 0r 1 f\\ L2 _ s(Rn) < |||||/|U 2( ^). (1-8) 

But since (— A — k 2 — iO)^ 1 is the integral operator of convolution type we can 
conclude that for fixed k it maps continuously H$ 1 (R n ) to H i _ s (R n ), where 
Hg l (R n ) denotes the dual of the Sobolev space H^ s (R n ). 
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We rewrite (1.6) as the integral equation 

u.sc = u + L k (u sc ), u = L k (u ), 
where the integral operator L k is defined as 

L k f(x) := J G+(\x - y\) (iV(W(y)f) + iW{y)Vf - q(y)f) dy. (1.9) 

The main result of present article is Theorem 2.1 which provides the exis- 
tence of the scattering solutions for the magnetic Schrodinger operator. This 
theorem is proved in Section 2. Based on the main result we justify in Section 
3 the direct Born approximation for such operators. 

2 Existence of the scattering solutions 

We are preceding a proof of the main result by the following lemmas. 

Lemma 2.1. Suppose that conditions (1.2) and (1.3) are fulfilled. Then 
there is 5 > \ such that ii e H\ 5o (R n ) and the integral operator L k maps 
H-s a (R n ) into itself. 

Proof. Conditions for p and S from (1.2) and (1.3) imply that there is 5 > \ 
such that 

L p s (R n ) C Ll(R n ). 

It is therefore true that under the conditions (1.2) and (1.3) functions V, W, 
VW and \W\ 2 belong to Lj o (R n ) with the same S - Since u is a bounded 
and smooth function we may conclude (using Agmon's result (1.8)) that 
uq belongs to H\ So (R n ). It can be mentioned here that we have no longer 
uniform estimates in k as in (1.8). In order to prove that L k maps H\ 5o (R n ) 

into itself we note that if / belongs to H\ Sa (R n ) then (1 + \x\ 2 )~~z f belongs 
to usual Sobolev space i^ 1 (i? n ). Using now Sobolev imbedding theorem we 
conclude that 

In 

feL^(R n ), n = 3, feLU (R n ), s < oo, n = 2. 

Then the conditions (1.2) and (1.3) and Holder inequality allow us easily con- 
clude that qf and {VW)f belong L 2 5o (R n ). Since we have W(x) G Lf(R n ) 
the function WV f will belong to Lj o (R n ) too. The final step is the applica- 
tion of Agmon's result (1.8). 
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We may prove a little bit more about this operator L k . 
Lemma 2.2. Let us assume that W G Lf(R n ), 

VW(x) G L p s (R n ), n = 3, 3 < p < oo; n = 2, 2 < p < oo, (2.1) 
w/jere 5 > ^ - -, and 

V(x) G L p oc (R n ), n = 3, 3<p<oo; n = 2, 2 < p < oo, (2.2) 

— * 

and i/ia£ and /love special behavior at the infinity such that 

\V(x)\, \W(x)\, \VW(x)\<-^-, \x\^oo, (2.3) 

where fi > 2 for n = 2, 3. Then the operator L k is compact in H ) i S() (R n ) for 
some S > \. 

— * — * 

Proof. Let us choose R > large enough and represent V, W and VW as 

v = v 1 + v 2 , w = m + W 2 , VW = V#i + vw 2 , 

where the supports of the functions Vi, W\ and Vlfi are included in the ball 
Br = {x G R n : \x\ < R}, but supports of the functions V 2 , W 2 and VW 2 
are included in the set {x G R n : |x| > i?}. Without loss of generality we 
assume in addition that the functions V 2 , W 2 and VW 2 are continuous and 
satisfy the conditions (2.3) for all \x\ > R. 

Conditions (2.1)-(2.2) and the previous considerations imply that 

iV(W 1 (x)f(x)) + mWVfix) - qi(x)f(x) G L 2 (B R ) 

for any function / G H^_ So (R n ). But L 2 (Br) is compactly imbedded in 
H~ l (B R ) and therefore in H^{R n ) for S > \. It remains to mention now 
that due to Agmon's result (1.8) operator (—A — k 2 — iO)" 1 maps continu- 
ously the space Hj^(R n ) to the space H\ 5o (R n ). Thus, part of the operator 

L k that corresponds to V±, W\ and VW\ is compact in H\ So (R n ) for So > \. 
Since outside the ball Br these functions V 2 , W 2 and VW 2 satisfy the con- 
ditions (2.3) we may firstly conclude using Sobolev imbedding theorem that 
for / G H\ So {R n ) we have 

iV(W 2 (x)f(x)) + iW 2 (x)Vf(x) - q 2 (x)f(x) e L 2 Q (R n ) 

for some 5q if and only if fi > 25q + 1. But under the conditions of Lemma 
2.2 this criterion is satisfied if So > \ is chosen appropriately. 
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Since the conditions (2.3) are satisfied we can find two sequences (f>j{x) G 
C™{R n \ B R ) and ^(x) (vector- valued) G C™{R n \ B R ) such that 

\\4>j ~ V 2 \\l™(r™\b r ) -»> 0, 11^- - W 2 \\l™{r™\b r ) -»> 0, 

|| V^- - VW 2 |U~(fl»\B«) 
as j — > oo for any 5' < /i. These approximation properties imply that 



\iV((W 2 - ^)/) + 2(^ 2 - Vi)V/ - (& - Qj)f\\ L j (Rn) 



< 



< c|||V(# 2 -^-)l + \W 2 -^\ + \q2-Qi\\\L% w\\f\\Hl lo (&) -> (2.4) 



as j — > oo, where Qj = \ipj\ 2 + <pj. Since we can choose 5 > | and /i from 
condition (2.3) such that 2S < \i then (2.4) means that part of the operator 
L k which corresponds to V 2 , W 2 and VW 2 is compact in H^ s (R n ) too. 

Lemma 2.3. Under the same assumptions as in Lemma 2.2 for any fixed 
k > and for any f G H\ 5o (R n ) with some 5q > | £/ie following asymptotical 
representation holds: 

ik\x\ l, n ~ 3 /• 

= en " / e"*^ (w(W(y)f) + iW{y)Vf - q(y)f) dy+ 



+ I TTa^T I ' M^°°> ( 2 - 5 ) 

X 2 



|x| 



Proof. In view of (1.9) one must study the behavior for \x\ — > oo of the 
function 



where i?i_ 2 denotes the Hankel function of first kind and of order IL ^-. In 



order to do that we take two cases, k\x — y\ > 1 and /c|x — y\ < 1. For the 
first case we use the behav 
(see, for example, |5j), i.e. 



first case we use the behavior of the Hankel function H„li for large argument 



#Sa(*) = c n -^= + O ( — ) , z^+oo. 

Z \Z2 
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So that (k > is fixed) we have that, 



e ik\x-y\ k ^ I l 

G+(k\x-y\) = c n - ^r + O rsr , |x| ^ +00. 

F-y| 2 \F-y| 
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Consider two subcases, < |x| a and \y\ > \x\ a , where < a < \ is a 
parameter. In the first case we have (since a < |) 



n— 1 . . n— 1 



= |x| (1 + 0(|x| a 1 )), \x\ ->■ +00. 
That is why we have for |y| < |a;| a with < a < | that 

as |x| — > +00. Substituting this asymptotic to the integral 

Si ■= J G+(\x - y\) (iV(W(y)f) + iW{y)Vf - q{y)f) dy 

k\x— y\>l 

gives that 

Si= J G+(\x - y\) (iV(W(y)f) + iW{y)Vf - q{y)f) dy+ 

k\x-y\>l,\y\<\x\°- 

+ f Gt(\x-y\) (iV(W(y)f) + iW{y)Vf - q{y)f) dy = 

k\x-y\>l,\y\>\x\ a 
^ifc I x I J^, 



= c n e |x| ^ J e ' tk(d '' y) (iV(W(y)f) + iW(y)Vf - q{y)f) dy+ 

k\x-y\>l,\y\<\x\ a 

+ J G+(\x-y\) (iV(W(y)f) + iW{y)Vf - q{y)f) dy+ 

k\x-y\>l,\y\>\x\° 

+ J 0(\x\^- 2a ) (iV(W(y)f) + iW{y)Vf - q{y)f) dy. (2.6) 

fc|x-j/|>l,|j,|<|x|° 

The conditions (2.1)-(2.3) allow us easily conclude that for any / e H\ 5o (R n ) 
the integrand iV{W{y)f) + iW (y)V / — q{y) f belongs to L l {R n ). Hence, the 
last term in the latter sum is o^x^ 11 ^) since < a < |. Denoting the second 



term in the latter sum (2.6) by I we may estimate it (using conditions (2.3)) 
as follows: 



\I\ < c 

k\x-y\>l,\y\>\x\ a 

c 

< 



\iV(W(y)f)+iW(y)Vf-g(y)f\ 
\x — y\ 2 



dy< 



\X\ 2 



+C 



iV(W(y)f) + iW(y)Vf - q(y)f\ dy+ 
\VW(y)\\f\ + \W(y)\\Vf\ + \q(y)\\f\ 



|x|«<| W |<Jfl 



\x — y\ 2 



dy< 



\X\ 2 



I/I + IV/I . 

,n^i, , dy. 

\x — y\ 2 | 1 m 



Since / e H_ So (R n ) the latter inequality implies 



l/l < o 



C 



I I n ~ ^ I I I n ~~ ^ 
X 2 / X 2 



/ 



y | n— 1 1 y 1 2/i— (n— 1) — 25o 



/ 



(2.7) 

Since /i > then 5 > | can be chosen here such that the last integral 
in (2.7) might be considered as the convolution of "weak singularities" and 
therefore we have 

' 1 ' (2, 



1 = 



\X\ 2 



The first case k\x — y\ > 1 is thus completely investigated. 

In order to consider the second case k\x — y\ < 1 we use the behavior of 
the Hankel function H^} 2 for small argument (see [5]) 



cz 2 ; n > 2, 
c(l+log(z)), ra = 2, 



z -> 0. 



Let us consider first n = 3. Then using this asymptotic and taking into 
account that for fixed k > and large |x| it can be assumed that \y\ > 



N 

2 ' 



one can estimate the integral 

S 2 := J G+(\x - y\) (iV(W(y)f) + iW(y)Vf - q{y)f) dy 

k\x—y\<l 
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as 



Sol < C 



J \Gi(k\x-y\)\ 2 dy 

\c\x-y\<l 



\ 



( 



i 

\~ 2 



|v/| 2 + |/p 



\y 



2fi 



dy 



< 



< 



C 



( 



\x 



H-5 



\ 



\X 



y\ 2 dy 



H 1 



s 



(R n ) 



J 



(2.9) 



\t\x-y\<l 

for \x\ — > +00, if So is chosen such that \ < So < 1. In the same way we have 
that 

/ 1 \ 

\x\ — > +00 



\x\ 2 



in the two dimensional case due to the behavior of the Hankel function 
for small argument. Combining (2.6), (2.8) and (2.9) we get (2.5). 

Remark 2.1. The proof of the last lemma shows that the function L k f(x) 
is continuous for all x such that \x\ > R, where R is large enough. 

These lemmas allow us to obtain the main result of this work. 

Theorem 2.1. Under the same assumptions as in Lemma 2.2 and for any 

k 7^ the integral equation (1.6) has a unique scattering solution from the 
space H\ So (R n ) for some So > \. 

Proof. Since the operator L k is compact in the space H^_ So (R n ) we can apply 
the Riesz theory in this Hilbert space. Based on this methodology we will 
prove that the homogeneous equation f—L k f = has only the trivial solution 
in the space H\ 5o (R n ) (i.e. the operator / — L k is injective). But this will 
imply that the operator is also surjective and the inverse (J— Lfc) 1 is bounded 
in H\ So (R n ). This condition is equivalent to the claim that the equation (1.6) 
has a unique solution from the space H\ 5 (R n ). 

It is possible to check that any / G H\ 5o (R n ) which satisfies the homo- 
geneous equation / — L k f = belongs to Hf oc (R n ) and satisfies also the 
equation 

Hf = k 2 f 

and Sommerfeld radiation condition (1.4). These facts imply that 
0= J v(j(V + iW)f-f(V-iW)j)dx. 

\x\<R 
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The divergence theorem then gives that 

0= f (JdJ-fdJ)da(x) + 2i J \f\ 2 Wuda(x), 

\x\=R \x\=R 

where v denotes the normal vector at the boundary of the ball B R . The 
radiation conditions (1.4) for the functions / and / allow us to conclude that 

= 2i y \f\ 2 Wvdo{x) + 2ik J \ffda(x) + o(^^J J fda(x). 

\x\=R \x\=R \x\=R 

Now the assumption that the function / e H\ SQ (R n ) satisfies the homoge- 
neous equation / — L k f = implies by the equation (2.5) that 

This behavior and the condition (2.3) imply that 

J \f\ 2 da{x) = o(l), R^oo. 

\x\=R 

This fact and Lemma 2.3 (see (2.5)) allow us easily conclude that (it is enough 
to integrate in (2.5) with respect to x) 

J e-af*) (iV(W(y)f(y)) + iW(y)Vf(y) - q(y)f(v)) dy = 0. 

Thus, we have actually (since the function / satisfies the homogeneous equa- 
tion / — L k f = 0) that 

/ = (j^f) ' |x| ^ °°- (2 - 10) 

The Sommerfeld radiation condition (1.4) implies also that 

V/ = o|— U), khoo. (2.11) 

\x\ 2 
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Using these two facts we are going to prove that / = a.e. In order to do 
so we first prove that f(x) = for \x\ > Rq with R large enough. Indeed, 
it suffices to prove that, for r > i? , the radial function 

F(r) := J f(r9)(f)(9)de (2.12) 

gn-l 

is identically equal to zero for each eigenfunction cf> of the Laplace operator 
A5 on the unit sphere S n ~ l . The eigenfunctions satisfy the equations 

(A s + /i 2 )0 = 0, fj 2 = k(k + n-2), k = 0,1,2,..., 

where /i 2 are the eigenvalues of the Laplace operator (see, for example, 

In view of the formula for the Laplacian A on R n in polar coordinates 

d 2 n — ld 1 

it follows that F(r) satisfies the ordinary differential equation 
F"{r) + — -F'{r) + (k 2 - ^)F(r) = 

= -2i J W{re)Vf(r6)<f){6)d6- J V{r6)f{r9)4>{6)d6, (2.13) 

where V = iVW + \W\ 2 + V. Let us rewrite this linear equation in the form 

F"( r ) + ^—If'M + (k 2 - ^r)F(r) = $(r, F). 

Since V and W satisfy (2.3) then it is not difficult to check that function 
$(r) from the right-hand side of equation (2.13) for r — > 00 behaves as 

$(r,F) = o(Jz), a>n. (2.14) 

It is well-known (see, for example, [5]) that the homogeneous equation cor- 
responding to (2.13) has two linearly independent solutions r~^~Hy\kr), 
j = 1, 2, where H^\z) and Hu\z) are the Hankel functions of order v, v 2 = 
/i 2 + (^^p) 2 , and of the first and the second kind, respectively. In view of 
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the asymptotic behavior of the Hankel functions (see |5]) it follows that the 
behavior of these two solutions is of the form 

r-^Hi>\kr) = -%e ±lkr + o(-^l), r oo, (2.15) 

f 2 \<f 2 / 

where + corresponds to j = 1 and — corresponds to j = 2. Next, we use 
Green's function (see, for example, [25J) for the Bessel equation with Dirichlet 
boundary conditions on the interval [r ,2r ], where r is large enough, 

Here u\ and u 2 are two linearly independent solutions of homogeneous equa- 
tion (2.13) which satisfy homogeneous Dirichlet boundary conditions at r 
and 2r , respectively, that is, 



u 2 lr 



r r 
1 



n-2 
2 



«i(r) = ( — ) (Hl l \kr )Hi 2 \kr) - H^\kr)H^\kr,)) , 



{H^\2kr,)H^{kr) - (kr) (2kr )) 



2r$r t 

and C = — Tg~ 1 'u 2 (ro)M / 1 (ro). Using this Green's function, equation (2.13) 
can be reduced to the following linear integral equation 

2r 
ro 

where K\ and K 2 are constants. This integral equation can be uniquely solved 
by iterations. Since (2.14) and (2.15) hold then using representation (2.16) 
we obtain for the integral part of the latter integral equation the following 
estimate: 

9(r,mZ,F)dZ = o(^ + -^^), a>n, (2.17) 



o 



as r — > +oo and r < r < 2r . Thus, the estimates (2.17) allow us to 
conclude that any solution F(r) of the non-homogeneous equation (2.13) has 
asymptotic (when r — > oo) that is the linear combination of two asymptotic 

(2.15). Since the hypothesis implies that F(r) = o ^-^-j, we deduce that 
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V(r) =0 for all r large enough. Thus, the same is true for f(x) (see definition 
(2.12)) for all \x\ > Rq with Rq large enough. 

We are in the position now to apply the unique continuation principle 
(UCP) in R n . Due to UCP for a second order elliptic differential operators 
with real coefficients (see Theorem 17.2.8 in [BJ) we may immediately con- 
clude that / = in R n . Thus, I — L/~ is injective and therefore, the integral 
equation (1.6) has a unique solution from the space H\ 5 {R n ) which is given 
by the formula 

Use = (I ~ L k y l u Q U = u + (I - L^LkUo, (2.18) 
where u is as in (1.5). Theorem 2.1 is completely proved. 

Corollary 2.1. If the conditions (2.1)-(2.3) are satisfied then the magnetic 
Schrddinger operator H has no positive eigenvalues. 

Proof. If A is a positive eigenvalue of H then 

Hu = Xu, ueL 2 (R n ), Hu G L 2 (R n ). 

It means that this u belongs to the domain of the Friedrichs self-adjoint 
extension of H and therefore u G W\{R n \ This fact allows us to conclude 
that this u satisfies the homogeneous equation 

u = L^u 

with an integral operator from (1.9). Thus, since W 2 1 (i? ra ) C H\ 6 (R n ) for 
5 > | (actually this imbedding holds for any 5 > 0) we may apply to this u 
the same proof as in Theorem 2.1 and conclude that actually u = 0. 

Using Agmon's results (1.8) the operator can be extended to the space 
L 2 _ s {R n ) as a uniformly bounded operator with respect to k > 1 such that 

IKIIz^) <C, k > 1, (2.19) 

where constant C depends only on the corresponding norms of V, W and 
VW. Based on this fact one can show that if the corresponding norms of V, 
W and VW are small enough then the operator norm of Lk as an operator 
from L 2 _ s (R n ) to itself is strictly less than 1. In that case the formula (2.18) 
can be rewritten as 

oo 

u = u + J2 L K u o)- (2- 20 ) 

i=i 

Thus, the scattering solution u can be obtained as the series of iterations of 
uq in the equation (1.5). 
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3 Scattering amplitude and direct backscat- 
tering Born approximation 

In this section we will consider the direct backscattering Born approximation 
for the magnetic Schrodinger operator H with conditions (2.1)-(2.3). The 
motivation to this problem is connected to the fact that the knowledge of the 
scattering amplitude with the backscattering data gives essential information 
about the unknown function V and W. 

Theorem 2.1 and Lemma 2.3 (see (2.5)) yield the following asymptotical 
representation for the scattering solutions u(x, k, 9) with fixed k > as \x\ — > 
+oo: 

ife | x |,"=3 / 1 \ 

u(x, k, 9) = e ik ^ + c n A(k, 9', 9) + o , 

\x\ 2 y\x\ 2 j 

where function A is called the scattering amplitude and defined by 

A(k, 9', 9) = J e - tk(e '' y) (iV{W{y)u) + iW(y)Vu - q(y)u) dy. (3.1) 

Substituting u = u + u sc into the equation (3.1) gives that 

A(k,9',9) = J e- lk{e '^ {iW(W(y)u )+iW(y)Wuo-q(y)u ) dy+ 

+ J e -W, v ) (^ V (w{y)u sc ) + iW{y)Vu sc -q{y)u s )} dy : = 

:=A B {k,9',9) + R{k,9',9). (3.2) 

The function A B is called the direct Born approximation. It can be checked 
(using integration by parts) that A B is actually equal to 

A B (k, 9\ 9) = -k(9 + 9')F{W){k(9 - 9')) - F(q)(k(9 - 9')), 

where F denotes usual n— dimensional Fourier transform as 

F(f)(0 = J /(■<■)< dr. 

The particular case 9' = —9 yields the direct backscattering Born approxi- 
mation 

A b B (k-9,9) = -F{q){2k9). (3.3) 
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Formulae (3.2) and (3.3) show that in the frame of the Born approxima- 
tion 

A(k, -9, 9) « -F(\W\ 2 + V)(2k9). 

But we want to write more terms in the Born series. For this purpose we 
calculate term R in the scattering amplitude. Using the series (2.20) and 
integration by parts we have that 

R(k,-e,6) = -i J e lk ^VW(y)L k u (y) dy+2k9 J e tk ^ y) W(y)L k u (y) dy- 

R n R n 

- J e ik ^q(y)L k u (y) dy + R 2 (k, -9, 9), 

R n 

oo 

where the term R 2 corresponds to the series L J k (uo) and equals to 

/oo 
e* k ^VW(y)J2 L lMy)dy+ 



3=2 



„ / oo \ „ oo 

+ 2ij e ik ^W(y)vlj2 L iMy)J dy - J e* k ^q(y)J2 L lMv) dy. 



R" \j~ 2 / 3-2 

(3.4) 

It will be shown that the term which correspond to R 2 in the definition (3.1) 
might be neglected because of the smallness of the operator norm L k in the 
space L 2 _ s (R n ). 

Since L k Uo(y) is equal to 

/ e *. fl ;,l s - i |,( i v# W -« W -,- w ),, 

R n 

then we obtain (after some simple calculations) the following representation: 
R(k, -9, 9) := R x (k, -9, 9) + R 2 (k, -9, 9) = 

= J J e tk ^ y+z) G+(\y- z\)VW{y)VW(z)dydz+ 

R n R n 

+Aik J J e ik{d ' y+z ^Gt{\y- z\)VW{y)9W{z) dydz- 

R n R n 
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-^k 2 J J e ik{d ' y+z) G+{\y - z\)9W(y)9W(z) dy dz+ 

R n R n 

+ 11 e lk ^y + ^G + k (\y - z\)q{y)q{z) dy dz + R 2 := 

R n R n 

:= h + I 2 + I 3 + I 4 + R 2 . (3.5) 

It can be mentioned here that this equality must be understood in the sense 
of tempered distributions. 

Using the facts F(G+)(v) = and W " = (^)~ n F(<P) * F{^) 

we can calculate the terms Ij,j = 1,2,3,4, more precisely as 

/l = (27f) J rfi-V-iO 

R n 



= Aik(27r)- n J 



F(VW)(k9 + rj)9F{W){k9 - rf) 



R n 



rf — k 2 — iO 



dr], 



I 3 = -4k 2 (27r)- n J 



9F(W){k9 + rj)0F{W){ke - rj) 
rf — k 2 — iO 



^ (2 ,r/ ^^)f«f-^ , (3,, 



R" 



Our next step is to neglect the term R 2 in (3.5) and justify this neglect. 
Indeed, using (1.7) and Agmon's results (1.8) the Z/^-norm of L k u can be 
estimated as (uniformly in \k\ > 1) 

\\LkUo\\ L 'i s (R") < c (\\VW\\ L 2 (Rn) + HW'llL^fln) + \\q\\ L 'j(R")) ■ 

The conditions (2.1)-(2.3) show that the right hand-side of the latter inequal- 
ity is finite. Thus, there is a constant Co depending only on the L^-norms of 
functions W , VW and q such that uniformly in |fc| > 1 

II^oIIl^^") < cq. (3.7) 

At the same time for any function / G H\ 5 (R n ) with some S > \ we can 
easily obtain 



\Lkf\\ L l s (R-) < C ( \\S7W\\ L V s{Rn) + ||W||L=o ( fl») + \\q\\ L p 2S (R«)) ||/||Hi 4 (fl»), 

(3.8) 
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where p is the same as in the conditions (2.1)-(2.2). We can rewrite (3.8) in 
the form of operator norm 

— * — * 

where constant C\ depends only on the norms of functions W, VW and q 
from (3.8). Hence, we may assume that these norms are chosen so small that 
Ci < 1. Now we extend the operator as an operator from L 2 _ 5 (R n ) to 
L 2 _ & (R n ) with the same norm estimate as in (3.9). This fact together with 
estimate (3.7) imply that 

oo 

nE L >oii^)^rrT- ( 3 - 10 ) 

Hence, the left hand-side of (3.10) can be made as small as we want if c\ (and, 
in addition, c ) are chosen small enough. This fact and duality arguments 
show that 

oo 

\\v(J2 l I u o)\\hzi^)<-t^i (s- 11 ) 

3=2 1 

The estimates (3.10) and (3.11) imply that one can have the term R 2 (k, —9, 9) 
as small as desired uniformly in \k\ > 1 and 9 G S 11 ^ 1 if the corresponding 
norms of functions W, VW and q (or the constants Cq and c±) are chosen 
small enough. Thus, the term R2(k, —9,9) can be neglected in the Born 
approximation. 

Summarizing our considerations (see (3.5)-(3.6) and (3.10)-(3. 11)) we may 
now obtain the following direct backscattering Born approximation (more 
precise than (3.3)) for the magnetic Schrodinger operator 

A(k, -9, 9) « -F(|H>| 2 + V)(2k9) + h + I 2 + I 3 + / 4 . (3.12) 

This formula gives us very good approximation for the backscattering am- 
plitude A. It is very important that for this approximation we need to have 

— * 

only the magnetic potential W and electric potential V, but we do not need 
(as we can see the formula (3.1)) to have the scattering solutions u(x,k,9) 
of the equation 

Hu(x) = k 2 u(x). 

This direct approximation (3.12) will be effectively used for the inverse 
backscattering Born approximation. Namely, due to formulas (3.6) we will 
be able to calculate precisely the quadratic term in the Born series that 
corresponds to the inverse backscattering approximation and to estimate its 
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smoothness. This smoothness result together with (3.12) will give us the 
solution of the inverse backscattering problem with respect to the recon- 
struction of the singularities and jumps of the unknowns. These problems 
will be investigated carefully in the futures publications. 
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